LETTERS 


these, 41.2% were characterized as GII. 3 
(mostly with a Gll.Pb pol), 26.5% as GII. 
Pg GII.l, and 17.6% as G1I.4 variants. 
From spring to fall 2012, the vari¬ 
ant New Orleans 2009 became the 
predominant GII.4 strain, and the vari¬ 
ant Sydney 2012 strain apparently 
disappeared. 

During November-December 
2012 and January 2013, ISGEV de¬ 
tected NoV infection in 90 (28.9%) of 
311 children hospitalized for gastroen¬ 
teritis. This finding is comparable to a 
prevalence of 25.2% in the same period 
(November-January) of the 2011-12 
winter season. A representative subset 
of 45 samples was randomly selected 
for sequencing, and 26 (74.3%) of 35 
fully typed strains were characterized 
as GII.4 Sydney 2012, which suggested 
that the new NoV variant had become 
the predominant strain in Italy. 

Our surveillance seem to mirror 
observations of a report from Denmark 
that documented the onset and circula¬ 
tion at low prevalence of the variant 
GII.4 Sydney 2012 at the beginning 
of 2012 with a marked increase in the 
prevalence only by the end of 2012 
( 10 ). Our surveillance detected the 
emergence of this variant in Italy at 
the end of 2011 and provided us with 
one of the earliest strains of the variant 
GII.4 Sydney 2012. This novel vari¬ 
ant has a common ancestor with the 
NoV GII.4 variants Apeldoorn 2008 
and New Orleans 2009 and has sev¬ 
eral amino acid changes on the main 
epitope in the capsid P2 domain ( 10 ). 

Sequence analysis of these early 
strains of the GII.4 variant Sydney 
2012 could help clarify the mecha¬ 
nisms driving its global emergence 
and spread. Continued surveillance 
for NoV infections through ISGEV 
and additional data on clinical and 
epidemiologic features will enable 
further assessment of the public health 
implications of the new variant GII.4 
Sydney 2012 in Italy. 
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Group C 

Betacoronavirus 
in Bat Guano 
Fertilizer, Thailand 

To the Editor: Bats play a criti¬ 
cal role in the transmission and origin 
of zoonotic diseases, primarily viral 
zoonoses associated with high case- 
fatality rates, including those caused 
by Nipah virus (NiV) and severe acute 
respiratory syndrome (SARS)-like 
coronavirus (CoV) infections ( 1 ). Re¬ 
cently, the World Health Organization 
(WHO) reported 44 confirmed cases 
of human infection with Middle East 
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respiratory syndrome CoV, resulting 
in 22 deaths. Full-genome and phy¬ 
logenetic analyses of these Middle 
East respiratory syndrome Co Vs 
have been published elsewhere (2). 
The identified viruses from 2 patients 
(previously referred to as England/ 
Qatar/2012 and EMC/2012) are ge¬ 
netically related and belong to group 
C betacoronavirus, which is most re¬ 
lated to Co Vs from Nycteris bats in 
Ghana and Pipistrellus bats in Europe 
(2,5). In addition, bat CoVs HKU4 
and HKU5 originated from Tylonyc- 
teris pachypus and Pipistrellus abra- 
mus bats, respectively, in the People’s 
Republic of China (4). Bats are also 
known to harbor and transmit nonvi- 
ral zoonotic pathogens, including the 
fungal pathogen Histoplasma capsu- 
latum, which causes histoplasmosis in 
humans (5). 

Bat guano is sold for use as a 
fertilizer in several countries, includ¬ 
ing Thailand, Indonesia, Mexico, 
Cuba, and Jamaica. The practice of 
collecting and harvesting bat guano 
may pose a considerable health risk 
because guano miners have a high 
level of contact and potential expo¬ 
sure to bat-borne pathogens. To as¬ 
sess pathogens in bat guano, we ex¬ 
amined bat guano from a cave in the 
Khao Chong Phran Non-hunting Area 
(KCP-NHA) in Ratchaburi Province, 
Thailand, where bat guano was sold 
as agricultural fertilizer, for the pres¬ 
ence of NiV, CoV, and H. capsulatum 
fungi. Bats from 14 species in 7 fami¬ 
lies have been found roosting within 
this area. Tadarida plicata bats are 
the most abundant species (2,500,000 
bats), and 3 other species of bats 
found at the site each had thousands of 
members: Taphozous melanopogon, 
Taphozous theobaldi, and Hipposide- 
ros larvatus. 

A random sample of dry bat gua¬ 
no, -100 g, was collected in a sterile 
plastic bag weekly from the main cave 
at KCP-NHA from September 2006 
through August 2007. The specimens 
were sent for analysis by express mail 


(at room temperature within 2-3 days) 
to the WHO Collaborating Centre for 
Research and Training in the Viral 
Zoonoses Laboratory at Chulalong- 
korn University. Samples were frozen 
immediately at -80°C until nucleic 
acids were extracted and PCR assays 
were run. A total of 52 collected bat 
guano specimens were examined in 
this study. 

Two aliquots of feces from each 
weekly specimen (104 samples total) 
were screened for CoV, NiV, and //. 
capsulatum by PCR. RNA was extract¬ 
ed from 10 mg of fecal pellet by using 
the QIAamp Viral RNA Mini Kit (QIA- 
GEN, Hilden, Germany). CoV RNA 
was detected by using nested reverse 
transcription PCR with the degenerated 
primers to amplify the RNA-dependent 
RNA polymerase (RdRp) gene (6). NiV 
RNA was detected by duplex nested re¬ 
verse transcription PCR (7). To detect 
H. capsulatum and other fungi, we 
extracted genomic DNA directly from 
bat guano by using the silica-guani¬ 
dine thiocyanate protocol, NucliSense 
Isolation Reagent (bioMerieux, Box- 
tel, the Netherlands), according to the 
manufacturer’s protocol. We tested 
for fungal ribosomal DNA (rDNA) in 
extracted total nucleic acid specimens 
by using the PCR protocol designed to 
amplify all rDNA from 4 major fungus 
phyla at the internal transcribed spacer 
1 and 2 regions (8). 

Four (3.8%) of 104 samples were 
positive for CoV. They were collect¬ 
ed on September 2, 2006 (KCP9), 
October 26, 2006 (KCP12), Novem¬ 
ber 14, 2006 (KCP15), and March 4, 
2007 (KCP31). Three of the 4 posi¬ 
tive CoV sequences (KCP9, KCP12, 
and KCP15) were identical at 152 nt 
of the RdRp region (ATCGTGCTAT- 
GCCTAATATGTGTAGGATTTTT- 
GCATCTCTCATATTAGCTC- 
GTAAACACAATACTTGTTG- 
TAGTGTTTCAGACCGCTtT- 
tATAGACTTGCaAACGAGTGT- 
GCGCAAGTCTTGAGTGAGTAT- 
GTGCTATGTGGTGGTGGCTAT) 
and phylogenetically clustered with 


the group C betacoronavirus (Fig¬ 
ure), with 76%, 80%, and 77% 
nt identity to bat CoV HKU4, bat 
CoV HKU5, and human CoV EMC 
and England ICoV, respectively. 
The other CoV sequence (KCP31: 
ATCGTGCACTTCCCAATAT- 
GATACGCATGATTTCCGCCAT- 
GATTTTGGGATCAAAGCATGT- 
TACTTGCTGTGACACATCT- 
GATAAGTATTACCGTCTTTGTA- 
ATGAGCTtGCACAAGTTTTGA- 
CAGAGGTTGTTTATTCTAATGG- 
TGGTTTC) showed 82% nt identity 
with bat CoV HKU8, an alphacoro- 
navirus. Although we recognize that 
longer sequences or full genomes 
may alter the topology of the phylog- 
eny slightly and give stronger branch 
support, we expect that the overall to¬ 
pology and placement of these CoVs 
would remain consistent. Samples 
from particular bat species could not 
be identified because bats of different 
species roost in this cave, and sam¬ 
ples were pooled during collection 
for bat guano fertilizer. The detection 
of CoVs in bat guano from the KCP- 
NHA cave in Ratchaburi was con¬ 
sistent with the previous finding of 
alphacoronavirus from Hipposideros 
armiger bats from the same province 
in 2007, but those researchers tested 
fresh bat feces (9). 

All bat guano samples screened 
by PCR were negative for NiV and 
Histoplasma spp. but were positive 
for group C betacoronavirus. The 
natural reservoir and complete geo¬ 
graphic distribution of this CoV are 
currently unknown. Although we did 
not isolate live virus from these sam¬ 
ples, the detection of nucleic acid and 
previous isolation of viruses from bat 
feces and urine (10) warrants some 
concern that guano miners might be 
exposed to bat pathogens in fresh ex¬ 
creta as well as in soil substances. We 
suggest that guano miners use pre¬ 
ventive measures of personal hygiene 
and improved barrier protection to 
reduce the possibility of exposure to 
zoonotic pathogens. 
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Figure. Phylogenetic tree of 3 coronaviruses (CoVs) isolated from bat guano collected 
in this study (KCP9, KCP12, and KCP15); 19 additional human and animal CoVsfrom 
the National Center for Biotechnology Information database are included. Construction 
of the tree was based on 152 nt of the RNA-dependent RNA polymerase gene region 
by maximum-likelihood method and GTR+I model with the 1,000 bootstrap resampling 
method implemented in MEGA5 (http://megasoftware.net/). Numbers on branches 
indicate percentages of bootstrap support from 1,000 replicates. The scale bar indicates 
the estimated 0.1 nt substitutions per site. HCoV-HKUIC, human CoV HKU1C 
(DQ415913); HCoV-HKUIA, human CoV HKU1A(DQ415903); HCoV-HKUIB, human 
CoV HKU1B (AY884001); MHV, murine hepatitis virus (NC001846); HCoV-OC43, 
human CoV OC43 (AY585229); PHEV, porcine hemagglutinating encephalomyelitis 
virus (DQ011855); BCoV, bovine CoV (AF391541); BatCoV-HKU9, Rousettus bat CoV 
HKU9 (NC009021); BatCoV/KW2E-F93/Nyc, Nycteris bat CoV (JX899383); BatCoV/ 
KW2E-F53/Nyc, Nycteris bat CoV (JX899384); BatCoV/KW2E-F82/Nyc, Nycteris bat 
CoV (JX899382); BatCoV HKU4, Tyionycteris bat CoV HKU4 (NC009019); BatCoV 
HKU5, Pipistrellus bat CoV HKU5 (NC_009020); HCoV-EMC, human betacoronavirus 
2c EMC/2012 (JX869059); England1_CoV, human betacoronavirus England 1 
(NC_019843); SARS_BatCoV HKU3, severe acute respiratory syndrome (SARS)- 
related Rhinolophus bat CoV HKU3 (DQ022305); SARSr-CoV-CFB, SARS-related 
Chinese ferret badger CoV (AY545919); SARSr-ciCoV, SARS-related palm civet CoV 
(AY304488); SARS_huSARS-CoV, SARS human CoV (NC_004718). An expanded 
version is online at wwwnc.cdc.gov/E I D/article/19/8/13-0119-F1.htm 
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